Metric of a tidally perturbed spinning black hole 



Nicolas Yunes-'^'Q and Jose A. Gonzalez^ 

'institute for Gravitational Physics and Geometry, 
Genter for Gravitational Wave Physics, Department of Physics, 
The Pennsylvania State University, University Park, PA 16802-6300 
^Theoretical Physics Institute, University of Jena, Max-Wien-Platz 1, 07743, Jena, Germany 
(Dated: Id: notes.tex.v 1.5 2005/07/29 13:40:21 yunes Exp ) 

We explicitly construct the metric of a Kerr black hole that is tidally perturbed by the external 
universe in the slow-motion approximation. This approximation assumes that the external universe 
changes slowly relative to the rotation rate of the hole, thus allowing the parameterization of the 
Newman-Penrose scalar tpo by time-dependent electric and magnetic tidal tensors. This approxima- 
tion, however, does not constrain how big the spin of the background hole can be and, in principle, 
the perturbed metric can model rapidly spinning holes. We first generate a potential by acting with 
a differential operator on i/jq. From this potential we arrive at the metric perturbation by use of the 
Chrzanowski procedure in the ingoing radiation gauge. We provide explicit analytic formulae for 
this metric perturbation in Kerr coordinates, where the perturbation is finite at the horizon. This 
perturbation is parametrized by the mass and Kerr spin parameter of the background hole together 
with the electric and magnetic tidal tensors that describe the time evolution of the perturbation 
produced by the external universe. In order to make the metric accurate far away from the hole, 
these tidal tensors should be determined by asymptotically matching this metric to another one 
valid far from the hole. The tidally perturbed metric constructed here could be useful in initial data 
constructions to describe the metric near the horizons of a binary system of spinning holes. This 
perturbed metric could also be used to construct waveforms and study the absorption of mass and 
angular momentum by a Kerr black hole when external processes generate gravitational radiation. 

PACS numbers: 04.25.Dm, 04.25.Nx, 04.30.Db, 95.30.Sf 



I. INTRODUCTION 

Gravitational wave observatories, such as LIGO and 
VIRGO, have the potential to study black holes in the 
strong field regime 0]. These black holes are expected 
to be immersed in a sea of gravitational perturbations 
that will alter the gravitational field of the background 
hole. Even though tidal perturbations are expected to be 
small relative to the background, they will be important 
in some astrophysical scenarios when attempting to pro- 
vide an accurate description of the non-linear dynamical 
orbital evolution of bodies around this background. The 
need for high accuracy in the description of the orbital 
evolution derives from the fact that gravitational wave 
observatories are extremely sensitive to the phase of the 
gravitational waves emitted by the system. Therefore, 
since this phase is directly related to the orbital evolu- 
tion, in some astrophysical scenarios it is necessary to 
take these tidal effects into consideration. 

The study of gravitational perturbations around Kerr 
black holes is important for several reasons. First, it 
is of astrophysical interest to study the flux of mass 
and angular momentum across a perturbed Kerr hori- 
zon h, 3, 4, Ij, which can be calculated through ma- 
nipulations of the tidally perturbed metric computed in 
this paper. Although this flux might be small for equal- 
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mass binaries, in extreme-mass ratio inspirals (EMRI) 
up to 5% of the total energy might be absorbed by the 
background hole. This absorption might slow down the 
orbital evolution increasing the duration of the gravita- 
tional wave signal Space-based detectors, such as 
LISA, will be able to observe and measure the gravita- 
tional waveforms of EMRIs, since they have particularly 
low noise in the low frequency band where such inspirals 
are common. Therefore, precise knowledge of the gravita- 
tional waveform including the tidal perturbations effects 
might be important in data analysis P|. Finally, the ex- 
plicit formulae of this paper might be useful to compute 
initial data near the horizons of a binary system of spin- 
ning holes. For example, Refs. 0, 23 make use of such 
explicit formulae for the non-spinning case to construct 
initial data via asymptotic matching. This data might 
be useful to the numerical relativity community because 
it derives from an approximate solution to the Einstein 
equation and, thus, we expect it to accurately describe 
the gravitational field of the system up to uncontrolled 
remainders. 

In this paper, we analytically construct explicit for- 
mulae for the metric of a tidally perturbed Kerr hole, 
where the perturbations of the external universe are as- 
sumed to vary slowly in a well-defined sense. Metric per- 
turbations for non-spinning holes have been studied in 
Refsjimilllillll using the Regge- Wheeler formal- 
ism [l5|. However, this method is difficult to implement 
for spinning holes because the metric will now depend 
on both radius and angle in a non-trivial way, render- 
ing the Einstein equations very difficult to solve. For 
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this reason, we use the Chrzanowski procedure |16l Il7| 
to construct the metric perturbation from the Newman- 
Penrose (NP) scalar ipQ. This procedure allows us to cal- 
culate the metric in the so-called ingoing radiation gauge 
(IRG), which is suitable to study gravitational perturba- 
tions near the outer horizon (r+), since there the metric 
is transverse and traceless. 

We will work in the slow-motion approximation, de- 
scribed in detail in Refs. |l|ll[ll|ll|ll|ll|ll|23, 
where we assume that the rate of change of the curva- 
ture of the external universe is small relative to the rota- 
tion rate of the background black hole, which in principle 
could be extremal. The external universe is completely 
arbitrary in that sense, as long as it respects the slow- 
motion approximation. For example, in the case where 
the external universe is given by a second black hole in a 
quasicircular orbit around the background hole, this ap- 
proximation is valid as long as their orbital separation is 
sufficiently large. In particular, this separation must be 
at least greater than the inner most stable circular orbit 
(ISCO)|2l|23, so that the binary is still in a quasicir- 
cular orbit. In that case, the curvature generated by the 
second hole would correspond to the external universe, 
which will change slowly as long as the orbital velocity is 
sufficiently small. In this sense, the slow-motion approx- 
imation will hold for astrophysically realistic binaries as 
long as they are sufficiently separated. 

Assuming this approximation to be valid, Poisson 2] 
has computed ■00 in the neighborhood of a spinning hole. 
First, the Weyl tensor of the spacetime is re-expressed 
in terms of the electric and magnetic tidal tensors of 
the external universe. Using the slow-motion approxima- 
tion, Poisson argues that these tensors will be spatially 
coordinate independent if evaluated at sufficiently large 
distances from the worldline of the hole. With this ten- 
sor, the asymptotic form of ipQ, denoted by V'Oj is com- 
puted far from the hole by projecting the Weyl tensor 
onto the Kinncrsly tetrad. This scalar will be a combi- 
nation of slowly varying functions of time, which will be 
parametrized via the electric and magnetic tidal tensors, 
and scalar functions of the spatial coordinates, which will 
be given by the tetrad. The asymptotic form of ipo now 
allows for the construction of an ansatz for tjjQ, which 
consists of its asymptotic form ipQ multiplied by a set 
of undetermined function of radius Rm{r). These func- 
tions must satisfy the asymptotic condition ^ 1 as 
r ^ r+, as well as the Teukolsky equation. This last 
condition is a differential constraint on Rm {r) which can 
be solved for explicitly, thus allowing for the full determi- 
nation of "i/jo- In this manner, the final expression for ?/'o 
is obtained and is now valid close to the horizon as well, 
in particular, approaching the perturbations generated 
by the external universe sufficiently far from the hole's 
worldline. 

Once ■00 has been calculated, we can apply the 
Chrzanowski procedure to compute the metric pertur- 
bation, still in the slow-motion approximation. This cal- 
culation contains two parts: the computation of a poten- 



tial (4') and the determination of the metric perturba- 
tion {hab) from ^ . In principle, one might think that it 
would be easier to try to compute hab directly from ■0o- 
Chrzanowski attempted this by applying a differen- 
tial operator onto ?/;o- However, Wald 23] discovered that 
doing this leads to a physically different gravitational 
perturbation from that represented by i/^o- Cohen and 
Kegeles 24] showed that by constructing a Hertz-like po- 
tential ^ from ■00 first and then applying Chrzanowski's 
differential operator to ^ instead leads to the real met- 
ric perturbation. This is the procedure we will follow to 
construct the metric perturbation in this paper. 

The construction of 4" requires the action of a fourth 
order differential operator on 0o, where here we follow 
Ori This potential is simplified by the use of the 

slow-motion approximation that allows us to neglect any 
time derivatives of the electric and magnetic tidal ten- 
sors. Once is calculated, we can apply Chrzanowski's 
differential operator to this potential ,2^ 2^]. In this 
manner, we compute the metric of a perturbed spinning 
hole from a tidal perturbation described by ^00 in terms 
of tidal tensors. These tensors are unknown functions 
of time that represent the external universe and which 
should be determined by asymptotically matching this 
metric to another approximation valid far from the holes 

The metric computed here, however, has a limited ap- 
plicability given by the validity of the slow-motion ap- 
proximation and the Chrzanowski procedure. The slow- 
motion approximation implies that we can neglect all 
derivatives of the tidal tensors. Furthermore, since we 
are working only to first non-vanishing order in this ap- 
proximation, it suffices to consider only the quadrupolar 
perturbation of the metric, since the monopolar and dipo- 
lar perturbations are identically zero. In perturbation 
theory, any I mode in the decomposition of the perturba- 
tion is one order larger than the I -t- 1 mode. Therefore, 
any higher modes or couplings of the quadrupole to other 
modes will be of higher order. The Chrzanowski proce- 
dure also possesses a limited region of validity, given by 
a region sufficiently close to the event horizon so that the 
spatial distance from the horizon to the radius of curva- 
ture of the external universe is small. This restriction 
is because the Chrzanowski procedure builds the metric 
as a linear perturbation of the background and neglects 
any non-linear interactions with the external universe. In 
particular, this restriction implies that this metric can- 
not provide valid information on the dynamics of the en- 
tire spacetime. However, if this metric is asymptotically 
matched to another approximation that is valid far from 
the background hole, then the combined global metric 
will describe the the 3-manifold accurately. 

We verify the validity of our calculations in several 
ways. First, we check that V'o indeed satisfies the Teukol- 
sky equation. After computing ^ , we also check that it 
satisfies the Teukolsky equation and the differential con- 
straint that relates to ^00. Finally, we check that the 
metric perturbation constructed with this potential sat- 
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isfies all of the Einstein equations to the given order. We 
further check that this perturbation is indeed transverse 
and traceless in the tetrad frame so that it is suitable for 
the study of gravitational perturbations near the horizon. 

This paper is divided as follows. Sec. ^1 describes the 
slow-motion approximation in detail, summarizes some 
relevant results from Ref. and establishes some nota- 
tion. Sec. mil computes from tpo while Sec. IIVI calcu- 
lates hab from 5'. Finally, Sec. Ivl presents some conclu- 
sions and points toward future work. In the appendix, 
we provide an explicit transformation to Kerr-Schild co- 
ordinates that might be more amenable to numerical im- 
plementation. 

In the remaining of the paper, we use geometrized units 
(G = 1, c = 1) and the symbol 0(a) stands for terms of 
order a, where a is dimensionless. Latin indices range 
from to 3, where is the time coordinate. The Ein- 
stein summation convention is assumed all throughout 
the paper, where repeated indices are to be summed over 
unless otherwise specified. Tetrad notation will be used, 
where indices with parenthesis refer to the tetrad and 
those without parenthesis to the components of the ten- 
sor. The relational symbol ~ stands for "asymptotic to" 
as defined in [2^, while the symbols <C and ^ are also 
to be understood in the asymptotic sense. In particular, 
note that if /(r) is valid for r <^ b, then this function is 
not valid as r — > 5. In this paper we have relied heavily 
on the use of symbolic manipulation software, such as 
MAPLE and MATHEMATICA. 



II. THE SLOW-MOTION APPROXIMATION 
AND THE NP SCALAR 

In this section we will describe the slow-motion ap- 
proximation in more detail and discuss the construction 
of "00 due to perturbations of the external universe. Both 
the slow-motion approximation and "00 have already been 
explained in detail and computed by Poisson in Ref. 0- 
Therefore, here we follow this reference and summarize 
the most relevant results for this paper while establishing 
some notation. 

Let us begin by discussing the slow-motion approxi- 
mation. Consider a non-spinning black hole of mass mi 
immersed in an external universe, with radius of curva- 
ture TZ. This external universe could be given by any 
object that lives in the exterior of the hole's horizon, 
such as a scalar field or another black hole. The slow- 
motion approximation requires that the external uni- 
verse's length scales be much larger than the hole's scales. 
In other words, for the non-spinning case we must have 
rni/TZ ^ 1, since these are the only scales available. 

For concreteness, let us assume that the external uni- 
verse can be described by another object of mass M^xt 
and that our hole is in a quasicircular orbit around it. 
Then, we have 

n rm + Mext ' V b ' ^ ' 



where V is the orbital velocity and b is the orbital sepa- 
ration. There are two ways of enforcing rrii/TZ <^ 1: the 
small-body approximation, where we let mi/Mext 1; 
and the slow- motion approximation, where V <ti 1. How- 
ever, in a future paper we might want to asymptotically 
match the metric perturbation computed in this paper 
to a post-Newtonian (PN) expansion j23| , which requires 
small velocities. For this reason, we will restrict our at- 
tention to the slow-motion approximation, which implies 
that we can only investigate systems that are sufficiently 
separated. The ISCO is not a well-defined concept for 
black hole binaries, but it has been estimated for non- 
spinning binary numerically |2lL |22| to be given approx- 
imately by 0^/500-^'^/ 10, where M and ijJisco are the 
total mass of the system and its angular velocity at the 
ISCO respectively. For spinning binaries the holes can 
get closer without plunging, where the value of the ISCO 
becomes a function of the spin parameter of the holes. 
Regardless of the type of binary system, the slow-motion 
approximation will hold as long as we consider systems 
that are separated by at least more than their ISCO so 
that they are still in a quasicircular orbit. 

The above considerations suffice for non-spinning holes 
because there is only one length scale associated with 
it, mi. However, for spinning ones we must also take 
into account the time scale associated with the intrinsic 
spin of the hole. The slow motion approximation does 
not constraint how large the spin of the background hole 
could be. However, in the standard theory, it is usually 
assumed that isolated holes will obey a < mi, where a 
is the rotation parameter of the hole and mi is its mass. 
If the previous inequality did not hold, then the hole 
would be tidally disrupted by centrifugal forces. When 
the background hole is surrounded by an accretion disk, 
however, some configuration might lead to a violation of 
the previous condition j2|| , but we will not consider those 
in this paper. Let us then define a dimensionless rotation 
parameter x — a /mi, which is now restricted to < x < 
1. The mass of the hole and this rotation parameter now 
define a new timescale, related to the rotation rate of the 
horizon and given by 

rH = ^=^-^{l + V^), (2) 

where fin is the angular velocity of the hole's horizon 
|29| . The slow-motion approximation then requires this 
time scale to be much smaller than the time scale asso- 
ciated with the change of the radius of curvature of the 
external universe. Text, *-e- Tn/Text ^ 1- This condition 
then becomes 

ml/7^<x, (3) 

since Text ~ T^- We take the above relation as the precise 
definition of the slow-motion approximation for spinning 
black holes. Following the reasoning that lead to Eq. 
we must then have Xi which, for a binary sys- 

tem, means that the orbital velocity of the system cannot 
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exceed the rotation rate of the background hole. Thus, 
the slow-motion approximation implies that, to this order 
and for spinning holes, we can also only consider systems 
that are sufficiently separated and where the black holes 
have relatively rapid spins. This reasoning does not im- 
ply that the Schwarzschild limit is incompatible with the 
slow motion approximation. Actually, Poisson 0] showed 
that when computing certain quantities, such as ipo, this 
limit can be recovered if we work to higher order in the 
slow-motion approximation. The precise radius of con- 
vergence of the slow-motion approximation to first order 
remains unknown, but an approximate measure of how 
small a separation the approximation can tolerate will be 
studied in a later section. 

Now that the slow-motion approximation has been ex- 
plained in detail, let us proceed with the construction of 
ipo, which is defined by 



(4) 



where Cabcd is the Weyl tensor of the spacetime and 
1°" and m° are the first and third tetrad vectors. This 
null vector m° is not to be confused with the m that 
we will introduce later in this section to denote the an- 
gular mode of the perturbation. Poisson works with 
the Kinnersly tetrad in advanced Eddington-Finkelstein 
(EF) coordinates, also known as Kerr coordinates, which 
are well-behaved at the outer horizon, given by r_|_ — 
nil + (to? — a^Y/'^. 

The calculation of ij^^ is based on making an ansatz 
guided by its asymptotic form far from the worldline of 
the hole but less than the radius of curvature of the exter- 
nal universe, i.e. r+ <§; r TZ. In this region, the Weyl 
tensor can be decomposed into electric £ab and magnetic 
Bab tidal fields, which are slowly varying functions of ad- 
vanced time V only. Thus, in this region, the only spa- 
tial coordinate dependence in -00 is given by the tetrad, 
namely 



^0 ~ -^z,„(w) 2Y2 



(5) 



where the tilde is to remind us that this quantity is the 
asymptotic form of tp^ and where 2Y^{6,(j)) are spin- 
weighted spherical harmonics, given by 
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±1/ 



±i4 



— sin6'(cos 9 ^ l)e 
i(l =F 2 cos 6* cos^ e)e^ 



(6) 



In Eq. lO, the quantities are complex combina- 

tions of the tidal fields given by Zm{v) = am{v) -\-i[3m{v), 



where 

ao{v) = £ii[v) + £22 (w), 

a±i[v) = fi3(w) T «f23(w), 

a-±2{v) = fii(u) - f22(w) T 2ifi2(w), 

Mv) = Bll{v)+B22{v), 

I3±i{v) = Bi3{v) =F iB2?.{v), 

I3±2{v) = Bii{v)-B22{v)T'^iBi2{v). (7) 

Note that in this region, ipQ is independent of radial co- 
ordinate. 

With Eq. at hand, Poisson makes an ansatz for the 
functional form of "001 namely 



^ z,„(u)i?,„(r)2K[ 



(8) 



where Rm{r) is an undetermined function of radius that 
must satisfy Rm (r) ^ 1 as r ^ r+ . If Eq. ||SJ is inserted 
into the Teukolsky equation, one obtains a differential 
equation for i?m(r). The angular part of the Teukolsky 
equation is automatically satisfied by the angular decom- 
position of ipo in spin-weighted spherical harmonics with 
eigenvalue E = 4 (Eq. 2.10 in Ref. ^^). The spatial part 
of the Teukolsky equation yields a differential constraint 
for Rm{r) namely, 

d 

x{l+x)— + [3{2x + 1) + 2imj] — 
ax-' ax 



2x + l 
x{l + x) 



R^{x) = 0, (9) 



where x is a rescaled version of the radial coordinate 
given by 



'r+ 



(10) 



and where the inner and outer horizons are given, respec- 
tively, by r± = m\± {ml — 0^)^/"^. We should note that 
r-(_ is an event horizon, while r_ is actually an apparent 
horizon. Solving this equation ^ one obtains 

i?m(r) ^ A^x-^{l+x)-'^F{-^,l,-l + 2imr,-x), (11) 

where is a normalization constant given by 

= -^TO7(1 + im7)(l + 4771272). (12) 
6 

In Eq. (fTT|l . the function F(a, 6, c; x) is the hypergeomet- 
ric function 



i^(a, 6, c; x) 



{a)n{b)n X 

n=0 



{C)r 



(13) 



where (a)„ is the Pochammer symbol defined as 

(a + 77- 1)! 



(a)„ a(a+ l)(a + 2) ... (a + 77- 1) 
(a)o = 1- 



(a-1)! 



(14) 
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For the present case, the series gets truncated at the 
fourth power and we obtain 



F(-4, 1, -1 + 2im7; -x) 
6 



1 



2imj — 1 
12 



(2zr7i7 — l)im'y {2im^ — l)im^{2im^ + 1) 
12 

(2zTO7 — l)im"f{2imj + l)(2im'j + 2) 



(15) 



where 7 is a constant given by 



(16) 



In this manner, Poisson calculates ipOj which encodes 
the gravitational perturbations of the external universe 
on a spinning hole in the slow-motion approximation. 
The full expression for iJjq is then given by 
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-x^ + 



12 



2im7 — 1 [2imj — l)imj {2im^ — l)im^ {2im^ + 1) 
12 



{2im"f — I)im7(2im7 + l)(2im'y + 2) 
I 



(17) 



where we have used the final abbreviation 



(18) 



without summing over repeated indices here, to group 
all terms that are spatially coordinate independent. We 
have checked that Eq. ((T7jl indeed satisfies the Teukol- 
sky equation [s^l for the s = 2 mode that corresponds to 
this scalar. Note that, in this final expression, summa- 
tion over the m = mode is removed because Poisson 2j 
has shown that it corresponds to the Schwarzschild limit 
and, thus, it contributes at a relative 0{mi/TZ) higher 
than all other modes. Also note that the final expression 
for ipQ only contains the quadrupolar I = 2 mode, once 
more because higher multipoles will be smaller by a rel- 
ative factor of rrii/TZ in the slow- motion approximation. 
For this reason, there are no mode couplings in the ipQ 
presented above. 

The time-evolution of the tidal perturbation will be 
exclusively governed by the electric and magnetic tidal 
fields. These tensors should be determined by asymp- 
totically matching the metric perturbation generated by 
this to another approximation valid far from the hole. 
However, in Refs. [3, E, 0, HO] it has been shown 
that when the external universe is given by another non- 
spinning black hole in a quasicircular orbit, these tensors 
scale approximately as 



m2 
If' 



(19) 



where TO2 is the mass of the other hole (the external 
universe) and b is the orbital separation (approximately 
equal to the radius of curvature of the external universe) . 
Note that the factor of in the denominator is neces- 
sary to make the tidal tensors dimensionally correct. Also 
note that in Eq. H19|l we have neglected the time depen- 
dence of the tidal fields, which generally is given by a 
trigonometric function, since we are interested in a spa- 
tial hypersurface of constant time. For the case where 



the other hole is spinning, Eq. Ijiyi) will contain correc- 
tions proportional to Xi but these terms will not change 
the overall scale of the tidal tensors. 

Throughout the rest of the paper we will generate plots 
of physical quantities, such as "00, 5* and hab- For the pur- 
pose of plotting, we will have to make two choices: one 
regarding the physical scenario that produces the pertur- 
bation of the external universe; and another regarding 
the parameters of the background black hole. As for the 
physical scenario, we will choose the external universe to 
be given by another orbiting black hole in a quasicircular 
orbit. This choice allows us to represent the tidal fields 
with the scaling given in Eq. H19|) . This scaling is not 
the exact functional form of the tidal fields and, there- 
fore, the plots generated will not be exact. However, this 
scaling will allow us to provide plots accurate enough to 
study the general features of the global structure of the 
quantities plotted, as well as some local features near 
the horizons. Regarding the parameter choice, we will 
assume toi = TO2, a — 0.99mi, and b = IOtoi, where 
M = mi + 7712 = 1 is the total mass. These choices are 
made in accordance with the slow-motion approximation, 
while making sure the system is not inside its ISCO. How- 
ever, the formulae presented in this paper should apply 
to other choices of physical scenarios and background pa- 
rameters as well, as long as these do not conflict with the 
slow motion approximation. 

Given the chosen physical scenario, the formulae in 
this paper should apply to other mass ratios and sepa- 
rations, as long as the orbital velocity does not become 
too large. For the background parameters chosen, the or- 
bital velocity is approximately V « 0.3, which indicates 
that, for fixed masses, we cannot reduce the orbital sep- 
aration by much more without breaking the slow-motion 
approximation. However, since we provide explicit an- 
alytic formulae for all relevant physical quantities, we 
can estimate their error by considering the uncontrolled 
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relation to obtain 




FIG. 1: Plot of the real part of the NP scalar i/jq along the 
a;-axis with the plotting parameters described in Sec. HTl 



remainders, i.e. the neglected terms in the approxima- 
tion. In particular, in a later section, we will see that 
the uncontrolled remainders are still much smaller than 
the perturbation itself even at 5 = lOM as long as we re- 
strict ourselves to field points sufficiently near the outer 
horizon of the background hole. 

In Fig^ we plot the real part of ipo with the plotting 
choices described earlier. Observe that as r becomes large 
the scalar asymptotes to a constant given by tpQ. Also ob- 
serve that the functional behavior of the scalar is drasti- 
cally different as r becomes small. In this figure, as well 
as in future figures, we have chosen to include an inset 
where we zoom to a region close to the horizons, so that 
we can observe its local and global behavior. For the or- 
bital parameter chosen, the inner and outer horizons are 
given by r_ w 0.43M and r+ « 0.57M. Observe from 
the inset in Fig.^that the scalar diverges at the horizons 
and as r ^ 0, which is due to the choice of tetrad. Fi- 
nally, observe that, except for where it diverges, the real 
part of tpo is of 0{V'^) in the entire 3-manifold. 



III. THE POTENTIAL 

In this section we will use ipo to construct the potential 
\l/ by acting some differential operators on the NP scalar. 
This potential must satisfy the vacuum Teukolsky equa- 
tion for the s = —2 mode together with the differential 
equation 



(20) 



where the overbar stands for complex conjugation and 
where the differential operator is given by D = I'^da 0, 

MM 

Ori 01 has shown that the above differential equations 
can be inverted with use of the Teukolsky-Starobinsky 



p 



where here A is given by 



A — r" ~ 2mir + a 



(21) 



(22) 



while p is a constant that for the time-independent case 
reduces to 

p = [l{l + l)-s^ + [l{l + 1) - s2 ^ |s| + 2]' . (23) 

In our case, since s and I refer to -00, s = 2 and ^ 2 so 
that this constant becomes p — 576. 

The differential operator 1?^ is given in spherical Brill- 
Lindquist (BL) coordinates {tBL,r,6,(j)BL) by 



(24) 



neglecting any time dependence, since time derivatives 
will only contribute at a higher order. In order to com- 
pute 5* in Kerr coordinates, we must transform the above 
differential operator. The transformation between Kerr 
and BL coordinates is given by 



( 2m\r 

dv = Mbl + avBL I — -. K 1 



dr = di^BL, 

dO = dOBL, 



obl -+- drBL 



(25) 



After transforming the differential operator I?^^ we ob- 
tain 



(26) 



because the drBL / d(j) term in the transformation cancels 
the (f) dependence. Note that r and 9 do not change 
in this transformation and, thus, A remains unchanged. 
Also note that is a scalar constructed from differential 
operators on ■f/^o and, since the latter is a scalar, ^' will 
also be gauge invariant. 

Before plugging in Eq. H17f) into Eq. H21|l to compute 
let us try to simplify these expressions. In the previous 
section, we defined the inner and outer horizons r_|_ and 
r_, as well as the new variable x. We can invert the 
definition of x so that it becomes a definition for r as a 
function of x and then insert this into Eq. (|22|l . We then 
obtain 



A = 47721(1 



(27) 



where we have defined r] = m\ — a?. It is clear now 
that when we combine the square of this expression with 
Eq. H17|l some cancellations will occur that will simplify 
all future calculations. 
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We are now ready to compute ^ , but first let us rewrite 
the function we want to differentiate, namely 

A^^o = ^™^(^) 2^2'"(^: ^bl), (28) 

where F{x) is shorthand for the aforementioned hyper- 
geometric function and where Cm is a new function of 
advanced time only given by 

Cm = -IGBmrf = ~16AmV^Zm{v). (29) 

Note that the angular dependence occurs in the spherical 
harmonics, while the only x dependence is in the hyper- 
geometric function. We can transform the operator 
to X space to obtain 



Re^P vsr 



nt — /9 — i f) 

U -Or - 2^1/2 Ox 




(30) 



FIG. 2: Plot of the potential ^ n along the x-axis with the 
plotting parameters described in Sec. UTI 



Applying all these simplifications, ^E* becomes 



* = ^ E ^"^ ^Y2"'{0A)^^dt[F{x)]. (31) 



We can now apply all derivatives to obtain 



(32) 



where we have used the shorthand f'^") which stands 
for the nth derivative of the complex conjugate of the 
hypergeometric function. This derivative is given by 



288 



(2im7 + l)im^{—2im'y + l){—2imj + 2) 



, (33) 



Note that we can reexpress the constant 7 in terms of 77 
as 7 = a/(277^/^). Interestingly this constant combines 
with Cm to return an overall constant that is purely real 
so that ^ is given by 

* = -^A2^r'"((?)e-™*z„(z;), (34) 



electric and magnetic tidal tensors. Decomposing ^ we 
obtain 



*7 



A 



= -"^ E ^'"("m(")cosm0-/3„(w)sinm(/)), 



A' 



^ 24 E ^™ ("m(w)sinTO0 + /3,„(t;)cos) 



(35) 



Note that the entire radial dependence is encoded in A^ , 
whereas the angular dependence is hidden in the spheri- 
cal harmonics. The time dependence occurs only in the 
the tidal fields that should be determined via asymptotic 
matching, as mentioned previously. This is the potential 
in Kerr coordinates associated with the tpo calculated in 
the previous section in the slow-motion approximation. 

In Fig.[21we plot "^fn with the plotting choices described 
in Sec. ^ Observe from the inset that the potential has 
nodes at both horizons. Also observe that the potential 
does not asymptote to a constant, but instead it grows 
quartically. This growth is due to the factor of A^ that 
dominates at large radius. Finally, note that the poten- 
tial is still of 0{V'^) for radii sufficiently close to the hole 
(roughly for r < 8m). 



where here stands for the I — 2 spherical harmonics 
with zero (j) dependence. We have checked that the poten- 
tial of Eq. H34|l indeed satisfies the definition of Eq. H20(l 
as well as the Teukolsky equation for the s = — 2 mode 
with angular eigenvalue E — 10. Note that "if has units 
of mass squared because the electric and magnetic tidal 
fields scale as the inverse of the mass squared. Further- 
more, note that eventhough -00 is singular at the horizon, 
^! is finite and actually vanishes there. 

Next, we proceed to decompose ^I* into real and imag- 
inary parts. The potential contains 2 complex terms, 
namely the (j) part of the spherical harmonics and the 



IV. THE METRIC PERTURBATION 

In this section we compute the metric perturbation by 
applying Chrzanowski's differential operator to the po- 
tential calculated in the previous section. The full metric 
of the spacetime is given by 

9ab = gab + iT-ab, (36) 

where gj^i^ is the background metric and hab is the per- 
turbation. Since ipQ and ^' were computed in Kerr coor- 
dinates, the background should also be in this coordinate 
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system. This background is given then by 



9oo 



B 
%1 

903 

B 
9l3 

922 



= - 1- 



= 1, 



--^(2asin2 e), 



-a sir? 6, 



2mi r 



{r' +a')sin'9+^^{a^smH), (37) 



where mi is the mass of the background black hole, a 
is its spin parameter, related to the angular momentum 
vector by 5* = ma, and where S = + cos^ 6. 

Let us now construct the metric perturbation. We will 
work with the form of the differential operator presented 
in Ref. namely 

hab = 2^{{-lak{S + a + 3f3-T){S + 4(3 + 3T) (38) 
-mamb{D - p){D + 3p) + ^(a'Wfc) [{D - 2ipi) 
((5 + 4/3 + 3r) + ((5 + 3/? - a - 7f - r) 
(D + 3p)]}vI/), 

where we have replaced p — p ~ ~2ipj. The metric 
perturbation constructed in this fashion will be in the 
ingoing radiation gauge (IRG), which is defined by the 
conditions 



hu — hin = hirn — him — ^mm — 0. 



(39) 



In Eq. (|38|l there are terms that depend on the dif- 
ferential operators D = I'^da and 6 — m°'da, which in 
turn depends on the tetrad. To be consistent, we will 
continue to work with the Kinnersly tetrad in Kerr coor- 
dinates given by 



,1,0- 



2a 



'(2) 



'(3) 



'(4) 



(0,- — ,0,0 



-\/2(r -|- ia cos 9) 



la sin 0,0, 1, 



smf 



(40) 



where the overbar stands for complex conjugation. Note 
that the m" vector is the same as the one in BL co- 
ordinates, but Z° and are different. The differential 
operators associated with this tetrad in Kerr coordinates 
are 



D = 
6 = 



Or + -^04, 
1 



\/2{r + ia cos 9) 



sin 9 



(41) 



where once more we neglect the time derivatives by use 
of the slow motion approximation. The covariant form 



of the tetrad in these coordinates is given by 
21] 

'a" 



L 



1. 



, 0, —a sin^ 9^ , 



(42) 



1 



= ^(A,0,0,- 



a sm t 

V21: 
1 

\/2S 



a cos 9, 0, 



-aA sin 



COS( 



[ar sin 9,0, aT, cos 9, — r sin 9{r^ 



-a')] 



where the superscript / and R stand for the imaginary 
and real parts respectively. One can show that if this 
tetrad is used we can recover the background metric 
with the formula = 2Z(-^nj-) + 2m(^g^ifiby 

Eq. contains terms that depend on the spin coef- 
ficients of the background [s^- These coefficients, also 
called Ricci rotation coefficients for the case where the 
tetrad is non-null, are simply contraction of the tetrad 
with its derivatives. In the tetrad formalism, these quan- 
tities can also be related to the Riemann tensor. Let us 
decompose the spin coefficients into real and imaginary 
parts, i.e. 



PR = 


r 

PI 


= -—cos9, 






Pr - 


cot 9, 

4 S 




a 

-coti 


?cos6'. 


TTfl = 


-v2— rrsm 


9 cos 9, 






TT/ = 


\/2 a / n 

-s^ - 


a^ cos^ 9) sin 9 




V2 
2 


a_R, = 


-v2— rsm 


9cos9~^'- 
4 S 


cot 9, 






V2 a . 
sm 9 

2 Y? ^ 


- a^ cos^ 9) 


V2 

4 


-cote 



— smf 



(43) 

These spin coefficients are the same as those obtained 
with the Kinnersly tetrad in BL coordinates. This in- 
variance is due to the fact that the spin coefficients are 
only tetrad dependent but still gauge invariant. 

We will split Eq. (|38|l into 4 terms in order to make 
calculations more tractable. The split is as follows: the 
term proportional to lah will be denoted term A; the 
term proportional to marine will be referred to as term 
B; the first half of the term proportional to l(^a'^b) will 
be denoted term C; and the remaining of this term will 
be referred to as term D. In this manner we have 

hab = hUfc((5 + a + 3/3-r)((5 + 4/3 + 3r)]^', 

/if, = [~mamb{D~ p){D + ip)]^, 

h'^b = [liamb){D~2ipj){5 + A(3 + 3T)]^, 

h^b = [l(aTnb){5 + i(i~a-^-r){D + ip)\^X^^) 



and the full metric perturbation is given by 



mh^b 



h^b)- 



(45) 
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With this split we can now proceed to simplify each erate with the differential operators. Expanding all terms 
expression more easily. One such simplification is to op- we obtain 



J 



hi, = 

hi = 

h^a, = 

hi = 



-IJb [5^* + * (4(5/3 + 3St) + (7/3 + 2r + a) 5* + (a + 3/3 - r) (4/3 + 3r) , 

l^amb) [DS^ + mDf3 + S'i'DT + (4/3 + 3r) - 2ipi (4/3 + 3t) * - 2ipiS-9] , 
li^arrib) [SD^ + 3^5p + 3p6^ + (3/3 - tt - r - a) + (3/3 - tt - r - a) 3p*] , 

I 



(46) 



These expressions give the metric perturbation in terms 
of the action of the differential operators on the spin co- 
efficients and the potential. 



A. Action of the Differential Operators 

In order to provide explicit formulae for the metric per- 
turbation we must investigate how the differential opera- 
tors act on the spin coefficients and on the potential. Let 
us first concentrate on the action of the differential oper- 
ators on the spin coefficients. After taking the necessary 
derivatives and decomposing the result into imaginary 
and real parts, we obtain 



Dp = Dp,R + iDpj = 
Dt 



— cos^ ( 



S2 



+ i 



„ ra 

2^cos( 



DP = D,s,R + iD0,i = - 



8p = 



6(3 



V2ra 

— ^ cot 6* cos ( 
V2 , 



Sp,R + iSp,i = — ^ ^ sin ^ cos ^ 



+i 



ar 



V2 



2 S2 



■ sm( 



cos^ e-r'^ + Zr^a^ cos^ 6* - 3 cos'' Or^a^ 



+ i 



ra cos t 



4S3 sin^ e 
^ 3a^ cos* 61 - cos^ 6 + Sr^ 
4sh?^E3 



—ra cos 6 



4 sin^ 
1 



St = 6r,R + i5r,i = -^a^ cos^ 9 



-r"^ +a^ cos^ 6 - 1a^ 



Y? 



1 -r"^ ^ a^ cos^ Q - 2a^ 
-ra cos 6» ^3 



(47) 



We now need to act with the differential operators on 



the potential itself. Wc can separate the dependence 
from these operators to obtain 



where 
D 



D^ 

{D5)^ 
((5D)4' 



Or — 2^— — 
A 



((5^)™* 



(48) 



ima 



2r?i/2 ^ 2?? x{l + x)' 

Sm = So (de + -r^) , 
V sm ft' / 

60 = So,r = -V26l (49) 

We can also compute the square of these operators acting 
on A^y™, where Y™ stands for the spherical harmonics 
with no (f) dependence. Doing so we obtain 



Z)^[A2] = 8(r - mif Am^a^ 
— Vlimair — mi), 

iasin 

— I- . 

S 



+ 
+ 
+ - 



V2E 



sin 2^ 1+2(52 



smt 



60m ( ia sin0 2 



+TTa sin26' 
sm^ V \/2S S 



sm 



■ (to — cos ( 



i6D)rn[A'Y 



5o (r^ 



TO 



-y™ [4A(r-TOi) 
smy / 

— 2iTOaA] , 
{DSUIA^Y""] = -V2SlA^ (r^ + 

+idDU[A^Y"^], (50) 



where the commas stand for partial differentiation. 

Wc now have all the ingredients to compute the action 
of the differential operators on the potential. Doing so 
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we obtain 

(-Dm)'* 



12 
24 



J2 >""^me-™'^[2(r-mi)-zma], 



smf 



7n\ —im(f) 



^ Y^e-'^^zra [2(r - mi)^ + A 
—rn?a? — 3ima{r — toi)] , 
~ 24 



sin 26* 





m sin 




V2I] 


TO \ 




+ 


(5oTO 


sine j 




sin0 



msin6' Sq ^ 



—a^ sin 261 



{6DU^ = 



\/2E 
-cose)]^™} 

[25a{r — mi) — Soima] 
V2, 



(5nTO 



sm 



(to 



sin 6* 



10^ X' 



m'-'Q , 



(51) 



In order to complete the calculation, we need to pro- 
vide explicit formulae for the first and second derivatives 
of the spherical harmonics. These derivatives are given 

by 



Yf^ = -2cos^6i + l±cos6', 



= sine (4 cos 61^1), 



Y 



±1 

99 
±2 



Y 



±2 



.99 



— smO {cos 6 =F 1) , 

9 1 1 
cos d± - cos9 + - . 
2 2 



(52) 



Note that the spherical harmonics and all of its deriva- 
tives are purely real. 



are always e"™''^ and Zm- The third term varies depend- 
ing on the differential operator. Let us define the third 
term with superscripts as 



= -j^Yl [2{r - mi) - ima] 



-T( 

24 ^ V 



+ — r™) So 

[2{r - mi)' + a 
w?o? — 3ima{r — TOi)] , 

a' sin 26* 



6 



smf 



zasmti 
721] 

1 fl 



JqTO 



_ 



— - — + —a' sin 26* 
V2S S 

cos6i)]y'"}, 

^ V sine 

[2(5o(r — TOi) — JoJTOa] , 

^ V2^2 ■y 1^ 



n— (to 

2 /) V 



sm 



12 



24 



J- fl 



^0^ 



(53) 



In general, if we want to decompose the product of 3 
complex quantities a, b and c we will obtain 



B. Decomposition into real and imaginary parts 



We will conclude this section by explicitly taking the 
real part of the metric perturbation, so as to have ex- 
plicit formulae for the metric in terms of only the real 
and imaginary parts of the spin coefRcients, the poten- 
tial and the action of the differential operators on these 
quantities. 

Before decomposing Eq. however, we must de- 

compose the action of the differential operators on the 
potential, i.e. Eq. (|51|l . Let us first note that the action 
of any differential operator on the potential always con- 
tains the product of 3 complex terms, the first 2 of which 



{abc)B. = cr {aB.bR - a/&/) - c/ (a/;6/ + aibn) , 
{abc)i = cr {aiibi + aibji) + cj {anbn ~ aibj) . (54) 



Zm, their real 



Since we have identified a = e"™"^ and b 
and imaginary parts are aji = cosmcj), aj — ~-smm(j), 
i>R = am and &/ = —Pm- Finally, if we further decompose 
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c we obtain 



and where 



-R 



JDS) 



6 



12 ^ 



m/0 



Cr = - 



^ II 



m 
sin 6 



i^r"[2(r-mi) + A 



A 

12 



(DS) _ (SD) 



V2 
24 



(SD) V2 ^y-^ f m \ 



24 



{51 



a sin 



sin 20 2m 



sm( 



- cos 
' 24 



(^0% - Slj) ( 



TOa sin 26 



sin6lS sin 6*2 



(m 



cos ^^)^ + So 



sin 20 2m 



sin 6 



asinO 



m 



ma^ sin 20 

sin0S ^ sin 02 
ma 



(m 



V2S 



(55) 



.^/2E . V2E 
From these equations it is simple to reconstruct the real 
and imaginary parts of the action of the differential op- 
erators on the potential by combining Eq. H55|l and H54|l . 
For example, the real part of -Dm^ is then given by 



■ cos 0. 



(56) 



(A„«')i? 



c^R^ {aRbR - ajbi) - c^j^' {aRbj + aiha) 
^ y™ [— 2(r — mi) {am{v) costo0 



SD) 



12 



—Pm {v) sin m(f)) + ma (a™ (w) sin mcj) 
+I3m{v) cos met))]. (57) 



We are now finally ready to get a final expression 
for the metric perturbation by taking the real part of 
Eq. I^nil- Doing so we obtain 
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-4(5/3,/ - 3(5r,7 - 5T//3i? + 4a7/3i? - iaRTi - 4afl/3/) + ((J™^')/?, [1(3r + ur) + ((5™*)/ (-7/3/ - 2r/ + a/) 
+ [{5.nfmR] , 

+ (D„*)/ [im^mipR + 2(mf - mimi)pi] + a^-^ [(m>^ - m'^mf){p] - p^) - im^mipRpi 
+Dp,R{mimi - m^m^) + 2m^miDpj] + 64-/ [(m>^ - m^m^)pRpi + m^miipj - p^) 
1 



-AD 



ka (mf) {[(i?(5)„*]^ + 2(<5„*)/p/} - m^) {[{D5)„,-^]^ - 2(,5„^')/?,p/} + (i^™^')/?, [4/3/jm^ - (4/3/ + 3r/) 
(An*)/ [4/3/jm^) + (4/3/ + 3T/)m^] + {[3A,i? + (8/3/ + 6t/)p/ + 4^?/,,/^,] + (-3A,/ 
+ 8p//3/j)m^)} + */ [(-3 A,/ - 4Ajj + 8piPR)m§^ - [3A.fl + (8/3/ + 6t/)p/ + 4Z?0,/j]to^J ) 

^)R{~fn{)Pi + '>T^b)PR) ^ ■^{Sm'i')i{ml^pR + m^^pi) + (Ari*)fl 
(-a/? + 3/3/?, - tt/j ) + ( Ai^')/ (-3/3/;, + tt/? + a/;,) + m^' (-3/3/ - a/ + t/ - tt/ 
{-ttrpr + Sp^R - Q!/?,p/? - a/p/ + 3prPr - nipi + piTj - 3p//3/) + to^) {ctRpi - Spj - 3/3/ p/? 
+7r/?,P/ - vT/p/?, - 3p//3/? + TipR - a/p/j)] + 3*/ m^^{-aipR + ttrpi - Spj - vr/p/?, + r/p/j, - 3/3/p/j 
+aRPi - SpiPr) + ml^{-5p^R + aipi + orpr + ttrPr - SPrpr + ipiPi - p/r/ + nipi) 

I 



™6) i-^i -^i +Ti - 3/3/) 



(58) 



The full metric perturbation is then given by 



b J 



ab 



(59) 



This is the metric of a tidally perturbed Kerr black hole 
in Kerr coordinates. We can transform this metric to 
Kerr-Schild coordinates, but this is left to the Appendix. 
We have checked that this metric indeed satisfies the Ein- 
stein equations by linearizing the Ricci tensor and veri- 
fying that all components vanish to first order. We have 
further checked that the metric perturbation is trans- 
verse and traceless (/i° = and habl"' = 0) in the tetrad 
frame making it suitable to study gravitational pertur- 
bations near the horizon. Furthermore, we have checked 
that the conditions that define the IRG [Eq. H39|) ] are also 
satisfied. Another feature of this metric is that its deter- 
minant is zero, which renders it non-invertible. However, 
the full metric gab = gab ~^ ^ab is invertible and, thus, the 
calculation of the Einstein tensor is straightforward. 

The metric perturbation has now been expressed en- 
tirely in terms of quantities explicitly defined in this pa- 
per. These quantities are the real and imaginary parts 
of the spin coefhcients, the potential and the action of 
the differential operators on the spin coefficients and the 



potential. The spin coefficients where decomposed in 
Eq. (|43l) : the potential was decomposed in Eq. (|35|) : the 
action of the differential operators on the spin coefficients 
is given in Eq. (|47|) ; and the action of these operators on 
the potential is decomposed in Eq. H54() and H55|l . 

The metric perturbation possesses the general global 
features that it diverges as r ^ oo and it either diverges 
or converges to a finite value as r ^ 0. The behavior as 
r ^ oo is to be expected because the Chrzanowski pro- 
cedure seizes to be valid far from the hole. On the other 
hand, the behavior as r — > is a bit more surprising. In 
this region there are two different types of behavior: ei- 
ther the perturbation remains finite or it diverges. These 
different types of behavior depend on the component and 
axis we are investigating. On the one hand, there are 
some components that either are finite and of 0{V^) or 
vanish as r — > for all angles, such as hoi, hn, hi2, hi^ 
and /i22. On the other hand, there are other components 
that diverge along certain axis as r — > 0. For example, 
^00, /i03 and /i33 diverge along the x-axis, the y-axis and 
the x-y diagonal, while ho2 and /123 diverge along the 
y-axis and x-y diagonal. This divergence is due to the 
choice of tetrad, since the fourth Kinnersly tetrad vec- 
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FIG. 3: Plot of the 00 component of the metric perturba- 
tion along the x-axis (solid line) and along the y-z diago- 
nal (dashed line) with the plotting parameters described in 
Sec. El 



tor clearly diverges as r — > 0. Note, however, that since 
the divergences occur well inside the inner horizon they 
will be causally disconnected with all physical processes 
ocurring outside the outer horizon and, thus, these diver- 
gences are irrelevant to most physical applications. This 
divergent behavior could nonetheless be avoided if a dif- 
ferent tetrad, such as the Hawking-Hartle one, is used to 
compute the perturbation, but this will not be discussed 
here further. 

In order to illustrate this global behavior, we have plot- 
ted /loo in Fig. |31 with the plotting choices described in 
Sec. m along the x-axis and the y-z diagonal {9 = 7r/4 
and (p = '''■/2). Observe that the perturbation diverges 
as r — > oo and as r — > along the x-axis, but it remains 
finite as r — > along the y-z diagonal. Everywhere else, 
and in particular near the outer horizon, the perturbation 
is of O(F^), where V = (M/&)i/2 is the orbital veloc- 
ity. Finally, observe that the perturbation vanishes close 
but not really at either horizon, remaining finite through 
them. 

The divergence of the perturbation can be used to ap- 
proximately determine the region of validity of the ap- 
proximation. As we have discussed previously and seen 
in Fig. 131 the perturbation will be valid inside a shell 
centered at the background hole. The inner radius of 
this shell can be approximately determined by graphi- 
cally studying the radius where hab > 0{V'^), approx- 
imately given by Tinner ~ O.OIM. This value for the 
inner radius of the shell is only an order of magnitude 



estimate, but is suffices to see that 



The 



causal structure of the region r < r_ is extremely com- 
plex and many of its features are known to be unstable 
under small perturbations that destroy the symmetries 
of the spacetime However, note that the perturba- 
tion is well behaved for r ~ <^ M and in particular it 
is of the order predicted by the approximation [0(1^^)]. 



In any case, the region r < r_ is hidden inside the event 
horizon and most astrophysical applications will be con- 
cerned with regions of small spacelike separations from 
the outer horizon, and not the inner one. 

The outer radius of the shell can be estimated by 
studying the fractional error in the perturbation, which is 
determined by comparing the perturbation to the uncon- 
trolled remainders in the approximation. The approx- 
imate error bars in Fig. (21 are given by an estimate of 
these uncontrolled remainders, Shah, which are due to 
truncating the formal series solution at a finite order. In 
the present case, this truncation is done at 0{V^)^ where 
the tidal fields provide this scaling. The uncontrolled re- 
mainders then will be of 0(V^) and should come from 
time derivatives of the tidal tensors, i.e. Shab oc hab- The 
argument of the tidal tensors is ujv, where uj — V/b is 
the angular velocity and where v is the advanced time 
coordinate. Any time derivative will pull out a factor 
of UJ, which will in turn increase the order of that term. 
However, in order for the uncontrolled remainders to be 
dimensionally consistent {Shab must have the same di- 
mensions as hab), we need to multiply the time deriva- 
tives of the tidal tensors by r, so that 



Shab ~ rhab ~ r{V/b)hab. 



(60) 



This line of reasoning, however, only leads to an order 
of magnitude estimate of the uncontrolled remainder. In 
principle, any dimensionless scalar function could be mul- 
tiplying this estimate as long as it does not change the 
scaling. For the case of non-spinning holes, the metric 
perturbation has been computed to 0{V^) 0], which 
allows us to compare these terms to Eq. (|60|l . This com- 
parison suggests that the multiplicative scalar function is 
roughly unity and, thus, unnecessary so that Eq. (|60|l is 
indeed a good approximation to the scaling of the uncon- 
trolled remainders. Clearly, this estimate is not the ex- 
act error in the approximation, which can formally only 
be determined if we know the exact functional form of 
the next order term. This estimate, however, is a physi- 
cally well-motivated approximation for the uncontrolled 
remainders. 

The outer radius of the region of validity of this ap- 
proximation can be approximately determined by study- 
ing the behavior of these error bars. From Fig. Owe see 
that the error bars become considerable large (approx- 
imately 50% as big as the perturbation itself) roughly 
at r = 13M. However, for r < AM the estimated error 
bars are less than 10% relative to the perturbation. This 
seems to indicate that, even at b = lOM where the slow- 
motion approximation begins to become inaccurate, the 
approximate solution is still valid sufficiently close to the 
outer horizon. This estimation of the fractional errors is 
by no means a formal proof of the existence or size of the 
region of validity of the approximation. However, this 
estimation does provide a strong argument that the ap- 
proximation is indeed valid sufficiently close to the outer 
horizon. If 10% fractional error is tolerable, which would 
correspond to neglecting terms of 0{V^), then the outer 
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radius of the shell is approximately given by router ~ 4M. 
Clearly, as b is increased, the slow-motion approximation 
will become more accurate and, thus, the perturbation 
will be valid inside a bigger shell with larger outer radii. 

The results of this paper are clearly valid close to the 
outer horizon of the background hole, thus allowing the 
study of physical processes of interest to the relativity 
community. For instance, we can use the perturbation 
presented here to construct initial data for a binary sys- 
tem near either hole. In this case, we are not interested in 
the behavior of the perturbation inside the inner horizon 
because that region can be excised and does not belong 
to the computational domain. Furthermore, the calcu- 
lations of angular momentum and mass flux across the 
horizon can still be performed because they only depend 
on the behavior of the perturbation near or at the outer 
horizon. 



V. CONCLUSIONS 

We have computed a tidally perturbed metric for 
a spinning black hole in the slow-motion approxima- 
tion. This approximation allows us to parameterize 
the NP scalar ipo [Eq. H17|) ] in terms of the electric 
and magnetic tidal tensors of the external universe. 
With this scalar we can then construct a potential 
[Eq. (|34|l ] . by applying certain differential operators 
to it. From this potential, we can then apply the 
Chrzanowski procedure to construct a metric perturba- 
tion [Eqs. (|S5 t .(Pt.(|T7 |) . f5^ .(|S5|.(P l) . CT . 

The metric is naturally computed in the ingoing ra- 
diation gauge and in Kerr coordinates, which are suit- 
able to study perturbation near the horizon due to its 
horizon penetrating properties. This metric is given ex- 
plicitly in terms of scalar functions of the coordinates 
and is parametrized by the mass of the background hole, 
its Kerr spin parameter and electric and magnetic tidal 
tensors. The mass and the spin parameters of the back- 
ground hole do not have to be small relative to each other, 
so in principle the metric presented here is capable of 
representing tidally perturbed extremal Kerr holes. The 
tidal tensors describe the time-evolution of the pertur- 
bation produced by the external universe and, thus, are 
functions of time that should be determined by matching 
the metric to another approximation valid far from the 
holes. 

The slow motion approximation constrains what kind 
of perturbations are allowed, thus limiting the external 
universe that is producing them. In this approximation, 
the radius of curvature of the external universe must be 
changing sufhciently slowly relative to the scales of the 
background hole. One consequence of this restriction is 
that the tidal fields must be slowly-varying functions of 
time, thus allowing us to neglect their time derivatives. 
In this sense, this approximation leads to a quasi-static 
limit, where in this paper we have calculated the first 
non-vanishing deviations from staticity. Another conse- 



quence of this approximation is that we can parametrized 
the perturbation in terms of multipole moments, where 
here we have only considered the first non- vanishing one 
(the quadrupolar perturbation). In perturbation theory, 
the / + 1 mode will be one order smaller than the I mode, 
which allows us to neglect the octopole and any other 
higher modes, as well as any mode beating between the 
quadrupole and higher modes. 

Due to these restrictions, if we allow the external uni- 
verse to be given by another hole in a quasicircular orbit 
around the background (a binary system), we are lim- 
ited to those whose separation is sufficiently large. In this 
case, we must have sufficiently large orbital separations so 
that the Riemann curvature produced by the companion 
is large relative to the scales of the background hole. In 
other words, this approximation will break if we consider 
systems that are close to their innermost stable orbit and 
ready to plunge. We have seen, however, that for a sep- 
aration of 6 = lOM, the metric presented here is valid 
inside a shell given approximately by 0.01 A/ < r < AM . 

The reason why the region of validity of the approx- 
imation is a shell can be traced back to the choice of 
tetrad and to the limitations of the Chrzanowski pro- 
cedure. The perturbation cannot be valid too close to 
the background hole because in that region the Kinner- 
sly tetrad, from which the perturbation was constructed, 
is also divergent. This perturbation is also divergent as 
r — !■ oo because the Chrzanowski procedure builds the 
perturbation as a linear expansion of the metric. In other 
words, we cannot analyze the dynamics of the entire 
spacetime with this metric, since its validity is limited 
to field points sufficiently close to the outer horizon of 
the background hole. However, if this metric were to be 
asymptotically matched to another metric valid far from 
the background hole, then the resultant metric would be 
accurate on the entire 3-manifold up to uncontrolled re- 
mainders thus, capable of reproducing the dynamics of 
the entire spacetime. 

This metric might be useful to different areas in general 
relativity. On the one hand it might be useful in the 
construction of astrophysically realistic initial data for 
binary systems of spinning black holes. Its importance 
relies in that it could accurately represent the metric field 
in the neighborhood of a black hole tidally disrupted by 
a companion. These tides are analytic and arise as a 
true approximate solution to the Einstein equations in 
the slow-motion approximation. Initial data, however, 
requires a metric that is valid in the entire 3-manifold. 
Therefore, in order to use these results as initial data they 
will first have to be asymptotically matched 9, lOj to a 
FN expansion valid far from the background holes [27l| . 
In this manner, initial data could be constructed that 
satisfies the full set of the Einstein equation, including 
the constraints, to a high order of accuracy. 

Another use for the metric computed in this paper re- 
lates to the fiux of mass and angular momentum through 
a perturbed Kerr horizon. This flux will be important for 
EMRIs, where the effect of tidal perturbations could be 
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large enough to lead to large fluxes, which in turn could 
affect the gravitational wave simal emitted by the sys- 
tem Recent investigations |3| have used a curvature 
formalism to compute this flux directly from ipo- How- 
ever, there exists a metric formalism to obtain this flux 
directly from the metric itself. An interesting research 
direction would be to compute this flux and compare to 
the results obtained with the curvature formalism. 

Finally, the perturbed metric computed here can also 
be of use to the data analysis community to construct 
gravitational waveforms. EMRIs are particularly good 
candidates to be observed by LISA, but such observa- 
tions require extremely accurate formulae for the phas- 
ing of the gravitational waves due to the use of matched 
filtering. Recently, Ref. ^ studied how to use and imple- 
ment a quasi-Kerr metric (a perturbed Kerr metric in the 
limit of slow rotation of the background hole) to detect 
EMRIs with LISA. A similar study could be performed 
with the perturbed metric computed in this paper, which 
can also describe rapidly rotating black holes. 

Future work will concentrate on performing the neces- 
sary asymptotic matching to shape this metric into use- 
ful initial data for numerical relativity applications. The 
matching procedure will provide expressions for the tidal 
tensors in terms of PN quantities, as well as a coordinate 
transformation between Kerr coordinates and the coor- 
dinate system used in the PN approximation. In this 
manner, a piece-wise global solution can be computed, 
which will contain small discontinuities inside the match- 
ing region that could be eliminated by the introduction 
of transition function. Since these discontinuities will be 
small due to the matching, the transition functions will 
not alter the content of the data to the order of the ap- 
proximation used. After the matching is completed, we 
will have obtained an approximate analytic global metric 
that will contain the tidal fields of one hole on the other 
near the outer horizon of the former, where these fields 
come directly from solutions to the Einstein equations. 



Appendix 



In this appendix we provide explicit formula for the 
transformation from Kerr coordinates to Kerr-Schild co- 
ordinate. This transformation can be found in Refs. [sj 
and is given by 



X — sin 9 (r cos (f> — a sin ( 
y = sin 9 (r sin(j) + a cos ( 



z 



r cost 



(61) 



The inverse transformation is given by 



= arccos — , 
r 



— arctan 



ry — ax 
rx + ay 



(62) 



Other useful relations are 
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sin6' 



1/2 



ry — ax 



a'^){x'^ + y'^)] 
rx + ay 



1/2 ' 



[(r2 -Ka2)(x2 -hy2)] 



1/2- 



(63) 



Note that these transformation reduce to the usual trans- 
formation from spherical polar coordinates to Cartesian 
coordinates in the limit a — > 0. 



The Jacobian of the transformation, A"b = dx'^/dx^ 
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is given explicitly by 



1\. y 

z 

A^, 



w 
w 

R^ + a" 
w 

R^-a^ 



2r2 



^(1 + 



_(^2_^2)-l/2 

rx + ay 
(r^ + a2)(a;2 + y^) 



/ 2 2\~l/2 
/ 2 2\~l/2 



z2 / i?2^^2 



y-x 



w 

ry — ax 



+- 



R^-a^ 



w J 2r 
rx + ay 



a — r- 



rx + ay 
ry — ax 



rx + ay 

ry — ax 
y-x- 



(r2 + a2)(x2 + 2/2) |^\^ rx + ay 

R"^ — a'^\ y ry — ax 

H ^ + r - a 



w J 2r 

rx + ay 



(r2 + a2)(a;2 + y"^) 
R^ + a^\ z 



rx + ay 

ry — ax 



y-x 



rx + ay 



(1 



(1 



(1 



+- 



w J 2r 



Note that this Jacobian reduces to the standard Jaco- 
bian of the transformation between spherical polar and 
Cartesian coordinates in the limit a ^ 0. There is a more 
elegant way to write this Jacobian in tensor notation as 



A"" 
A^ 



= ^ ^5aix' + ^ [5iax' {R^ - a^) + 2a^zSaz] I . 



1 A'^ 

1/2 ^' 

r (r2 — a2) 



rx + ay 



(r2 + a2)(a;2 + y'^) 
ry — ax 



A\{y-'-^-^x] (65) 



+ r- 



rx + ay 



a 6ya 



rx + ay 
ry — ax 
rx + ay' 



r 4a 



With this Jacobian, the metric in Kerr-Schild coordi- 
nates is given by 



9ab = Qa'b'A'^ qA'' fe, 



(66) 



where here the primed indices refer to spherical coordi- 
(64) nates and the unprimed indices to Cartesian coordinates. 
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